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1 . INTRODUCTION 


There  is  a program  at  the  Harry  Diamond  Laboratories  to  begin  a more 
intense  testing  and  analysis  of  fluidic  component  characteristics. 
Because  of  the  magnitude  of  the  effort  required  to  obtain  precise  data 
on  fluidic  elements  ^ an  effort  has  been  made  to  standardize  and  automate 
as  much  of  this  procedure  as  possible.  This  effort  has  resulted  in  the 
construction  of  a fluidic  test  facility  that  includes  an  automated  data 
acquisition  system  and  computer  control  of  test  conditions.  A computer 
program  titled  "Fluid*'  has  been  written  to  perform  quasi-static 
performance  data  reduction  and,  specifically,  to  provide  the  information 
needed  to  make  up  the  digital  component  fact  sheet  adopted  by  the 
Government  Fluidics  Coordination  Group.  In  addition,  a standard 
interconnection  pattern  for  certain  in-house  components  has  been 
advanced. 

This  report  describes  CALFLO,  a program  written  for  the  WcUig  2200B 
programmable  calculator,  which  can  be  used  to  obtain  calibration 
information  on  Vol-0-Flo  meters.  The  input  to  the  prograim  is  the  data 
collected  with  the  HDL  flow  meter  calibration  equipment.  The  progrcim 
output  includes  a printout  of  the  data,  the  coefficients  of  the  second 
order,  least -squares  curve  fits  of  the  corrected  data  for  both  flow 
directions,  an  analysis  of  data  error,  and  plots  of  both  normalized  data 
and  calibration  curves. 


2.  DERIVATION  OF  EQUATIONS 

The  Vol-0-Flo  meter  is  a transducer  that  provides  a differential 
pressure  signal  (AP)  that  is  a function  of  the  volume  flow  passing 
through  its  midpoint.  The  governing  equation  for  this  transducer^  is 
given  by 


AP  = by  Q + CP  Q^ 
mm  mm 


(1) 


where 

b,c  are  constants, 

p is  kinematic  viscosity, 

Q is  volxime  flow, 

p is  density, 

m refers  to  conditions  at  the  midpoint  of  the  meter. 


I.  Steele  and  L.  B.  Orbacb,  Theory  and  Use  of  the  Vol-0-Flow 
Meter,  National  Instrument  Laboratories,  Inc.  (January  1963). 
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The  first  term  on  the  right-hand  side  of  equation  (1)  is  due  to 
Poiseuille  flow.  The  second  term  is  due  to  orifice  effects.  For 
Iconinar  flow,  the  first  term  is  much  larger  than  the  second. 

Assuming  steady  flow,  the  continuity  equation  for  the  mass  flow 
through  the  meter  and  the  calibration  system  is 

P Q = PQ  = M , (2) 

m m 

where  M is  mass  flow  rate. 

A third  equation  is  the  ideal  gas  law. 


where  T is  temperature  and  p is  pressure. 

When  a Vol-0-Flo  meter  is  calibrated,  Q , the  volume  flow  at  the 
meter  midpoint,  is  not  measured  directly.  Instead,  a known  volume  of 
gas  is  collected  at  a constant  flow  rate  over  a measured  period  of  time. 
However,  it  is  possible  to  calculate  with  equations  (2)  and  (3)  , if 
the  appropriate  pressures  and  temperatures  also  are  known.  This 
calculation  is  as  follows : 


where 

V is  the  known  volume, 
t is  the  measured  time, 

V refers  to  conditions  at  the  known  volume. 
Equation  (2)  gives 


pT/p  = constant 


(3) 


(4) 


(5) 
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Equation  (3)  gives 


P P T 
V _ V m 

P ” P T 

m m V 


(6) 


Assuming  that  T = T , 
^ m V 


then  combining  equations  (5)  and  (6)  gives 


(7) 


From  equation  (1)  , it  is  seen  that  for  a given  value  of  Q^,  Ap  is  a 
function  of  temperature  and  pressure,  since  \i  and  p are  temperature  and 
pressure  dependent.  Since  the  (CpQ^)  term  is  very  small,  the  changes 
in  p can  be  neglected  over  moderate  temperature  and  pressure  ranges. 
This  approach  allows  for  the  absorption  of  the  temperature  and  pressure 
corrections  of  y by  the  AP  term. 

To  correct  the  measured  value  of  AP  for  the  temperature  and  pressure 
variations  of  y,  the  following  derivation  is  used.  Since  AP  is  a linear 
function  of  y, 

AP  = AP  (Pg  + Ap)  = AP(Pg)  + Ap  , (8) 


where 

AP(ys  + Ay^  is  the  measured  value  of  AP  at  some  value  of  y = y^  + Ay, 

y is  the  value  of  y at  standard  conditions  of  temperature 
^ (20®C)  and  pressure  (1  bar) , 

Ay  is  the  difference  between  y and  the  actual  value  of  y. 

s 

From  this  discussion  and  equation  (1) , 


3(AP) 

3y 


= bQ, 


m 


(9) 
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Assuming  that  any  deviations  of  the  value  of  u from  a standard  value  can 
be  approximated  by  a linear  function  and  using  the  chain  rule,  one 
obtains 


(10) 


where 


AT  is  T - T , 
s in 

Ap  is  p - p , 
s m 

s refers  to  standard  conditions  (20*C,  1 bar) 

Combining  equations  (8)  to  (10) , one  obtains 

Ap  I 


AP  = AP  + 
s 


y9T  3p 


(11) 


where 


This  can  be  rearranged  to  give 


AP  = 
s 


AP 


, 1 1 3y  . 

14-  — *^AT  + — Ap 

3T  y^  3p 

s s 


(11a) 


_ , - 1 3y  _ , 1 3y  . _ 

The  values  of  — — AT  and  — — Ap  for  various  gasses  are  given  in 


y 3T 
s 


3p 


table  I. 
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TABLE  I.  SPECIFIC  GRAVITY  AND  VISCOSITY  AND  THEIR  CHANGES  WITH 
TEMPERATURE  AND  PRESSURE  FOR  VARIOUS  GASSES 


Gac 

Specific  gravity 
(Air  - 1) 

Ug  (N-»/ib*)  X lO’ 
at  T - 20*C,  p - 1 bar 

% Change  in  m 
for  change  in  T of  1*C 

m 

% Change  in  u 
for  change  in  p^  of  1 atm 

Air 

1.000 

161.5 

0.263 

0.13 

Argon 

0.5963 

225.2 

0.274 

0.14 

Carbon  dioxide 

1.5290 

146.5 

0.311 

0.25 

Helium 

0.13604 

195.4 

0.227 

0.011 

Hydrogen 

0.06952 

66.22 

0.230 

0.016 

Nitrogen 

0.96724 

175.4 

0.256 

0.13 

Oxygen 

1.10527 

202.7 

0.259 

0.13 

Equation  (1)  can  be  rewritten  as 

AP^  = BQ^  + CQ^  , (12) 

where 

B = by  r 
s 

C = cp^  , 


The  calibration  equation  now  becomes 


AP  = A + BQ  + CQ^ 
s mm 


(12a)l 


where  A is  the  zero  offset  of  a least-squares  curve 
This  can  be  written  also  as 

-B  + (b^  - 4Ci 


fit  to  the  data. 


2m  = 

m 


(b^  - 4c(a  - APJ)’’ 


2C 


(12b) 


flow  meter  is  not  at  standard 


of  APg  is  not  measured  at  the 


If  the  gas  at  the  midpoint  of  the 

conditions  (20®C,  1 bar)  ^ then  the  value 

signal  ports  of  the  flow  meter.  However,  equation  (11a)  can  be  used  to 

obtain  AP  from  the  measured  value, 
s 

The  value  of  APg  can  be  used  to  obtain  the  volume  flow  rate  at  the 

meter  midpoint  conditions,  Q^,  by  visual  use  of  the  plot  of  the 

calibration  curve  fit  or  by  use  of  the  value  of  AP  in  equation  (12b)  . 

s 
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The  volume  flow  rate  that  would  exist  at  any  arbitrary  conditions  of 
temperature  and  pressure  can  be  obtained  by  combining  equations  (5)  and 
(6)  : 


p T. 

®i  p . T ®m 
1 m 


.(13) 


where  i refers  to  any  arbitrary  conditions  of  T and  p. 

Equation  (12a)  is  used  to  obtain  the  least-squares  curve  fits  for 
the  corrected  calibration  data  on  the  flow  meter. 


3.  PROGRAM  INPUT  AND  OUTPUT 


This  program  was  written  in  BASIC  for  a Wang  2200B  calculator  with 
12k  core.  The  peripherals  used  include  a keyboard,  cathode  ray  tube 
(CRT)  display,  80-column  thermal  printer,  and  analog  flatbed  plotter 
10  X 15  in.  The  plotted  output  requires  graph  paper  10  by  15  in.  with 
the  plotter  zero  set  at  the  lower  left  corner  of  the  grid  and  the 
plotter  scale  set  8 in.  to  the  right  and  8 in.  above  that  point. 
Appendix  A lists  the  program. 

The  data  used  by  this  program  are  entered  from  the  keyboard.  The 
user  is  prompted  to  enter  the  data  by  the  CRT  screen.  The  units  in 
which  the  data  should  be  entered  (included  in  the  prompt  statements)  are 
meant  to  reflect  not  any  standard  \init  system,  but  rather  the  units  in 
which  the  instrumentation  used  to  calibrate  the  flow  meters  is  scaled 
(fig.  1). 

To  provide  a header  label  for  program  output,  the  program  asks  for 
specification  of  the  date,  meter  serial  number,  and  meter  range.  The 
meter  range  value  is  used  also  to  determine  the  scale  value  of  the  plot. 
Since  the  Vol-0-Flo  meter  is  usually  fairly  linear  up  to  twice  its 
nominal  flow  range,  useful  data  can  be  obtained  for  flow  levels  up  to 
twice  the  nominal  range.  Therefore,  the  plots  of  the  calibration  curve 
fits  are  automatically  scaled  to  twice  the  nominal  range  of  the  meter. 

The  calibration  data  are  entered  a point  at  a time  in  the  following 
order : 

(a)  Left  manometer  reading  (cmwc) 

(b)  Right  manometer  reading  (cmwc) 
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(c)  Back  pressure  (kPa) 

(d)  Volume  (liters) 

(e)  Time  (min) 

(f)  Temperature  (®C) 

(g)  Volume  pressure  (inwc) 

(h)  Barometric  pressure  (mm  Hg) 


VOL-O-FLO 


U TUBE 
MANOMETER 


Figure  1.  Schematic  of  HDL  flow  meter  calibration  equipment. 


Before  any  data  are  entered  from  the  keyboard,  a data  point  is 
automatically  generated  by  the  program.  This  point  corresponds  to  the 
condition  AP  = Q = 0.  As  the  values  for  each  complete  data  point  are 
entered,  they  are  printed  out  on  the  thermal  printer.  This  printing 
allows  the  user  to  check  his  input  data  before  program  calculations 
begin. 

If  an  incorrect  value  is  entered,  a correction  can  be  made  after  all 
the  data  for  the  current  flow  direction  are  entered.  The  program 
recognizes  that  the  last  datum  has  been  entered  when  the  value  for  the 
left  manometer  reading  is  entered  as  9.999.  At  that  point,  no 
additional  data  values  can  be  entered,  but  those  already  entered  can  be 
changed . 
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Changes  are  made  in  the  following  manner:  The  CRT  provides  a 

prompt,  asking  if  changes  are  to  be  made.  If  so,  it  prompts  again, 
asking  for  the  data  point  number,  which  was  printed  out  with  the  data  on 
the  printer.  (Its  value  is  one  higher  than  the  order  in  which  the  data 
were  entered,  since  the  program  generated  the  first  point 

for  AP  = Q = 0.)  The  program  then  asks  for  all  the  values  for  that  data 
point,  as  listed  above,  and  prints  them  out  for  subsequent  verification. 
The  CRT  again  provides  a prompt,  asking  if  changes  are  to  be  made.  This 
sequence  continues  until  all  changes  have  been  made. 

At  this  point,  the  progreun  calculates  a second-order 
least-squares  curve  fit  to  the  data.  It  then  prints  a data  analysis 
table  comparing  the  individual  data  point  to  that  curve  fit.  The  teible 
includes  the  following  parameters: 

(a)  AP 

s 

(b)  (calculated  from  data) 

(c)  0 (calculated  from  curve  fit) 


(d) 

Difference  between  0 from  data 

(b)  - (c)  ^ 

cind  Q 

from 

curve 

fit. 

(e) 

Percent  ratio  of  difference  between 
curve  fit  to  from  data,  [(d) /(b)j 

S' 

from 

100 

data 

and 

2,n 

m 

from 

(f) 

Percent  ratio  of  difference  between 

Q from  data 

and 

Q_ 

from 

curve  fit  to  full  range  flow  of  flow”'  meter,  [(d)/R]  100, 
where  R = range  of  flow  meter. 

The  program  next  prints  out  a value  for  the  sum  of  the  squares  of 
the  deviation,  based  on  (d)  . At  this  point,  the  user  can  look  at  the 
data  analysis  table  (column  f)  and  see  if  any  data  points  diverge 
significantly  from  the  curve  fit.  If  there  are  such  points,  the  user 
should  check  that  the  correct  data  were  entered.  An  entry  error  can  be 
corrected,  since  the  program  again  asks  for  data  corrections  at  this 
point.  Since  an  error  could  have  been  made  in  taking  the  data,  the 
program  allows  the  deletion  of  up  to  three  data  points,  once  all  desired 
data  corrections  have  been  made.  Once  these  corrections  and  deletions 
are  made,  a new  curve  fit  is  calculated,  and  a second  data  analysis 
table  is  printed.  The  program  continues  to  cycle  through  this 
correction,  deletion,  curve  fit,  and  data  analysis  routine,  until  the 
user  makes  no  corrections  or  deletions.  When  this  routine  is  completed, 
a table  is  printed  of  AP  versus  Q for  nominal  values  of  AP  and  values  of 
Q calculated  from  the  curve  fit.  The  program  also  provides  plots  for 
both  'AB'  and  *BA*  flow  directions  of  the  calibration  curve  fit  and  for 
the  corrected  data  on  which  it  is  based.  *Examples  of  program  output  are 
shown  in  figure  2 from  the  thermal  printer  (*AB*  direction  only)  and  in 
figure  3 from  the  plotter. 
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DATE  (DAY  MONTH  YR)  31/7/75 

METER  S/N  6014 

METER  RANGE  8.33  CCPS  (.5  LPM) 


FLOW  DIRECTION  IS  AB 


NO 

LEFT 

RGHT 

P2 

VOL 

TIME 

T2 

PO 

BARP 

DIFF 

1 

0.00 

0.00 

0.00 

0.0 

1.0000 

0.0 

0.0 

760.0 

0.00 

2 

3.29 

4.87 

1.50 

0.2 

1.3840 

22.1 

1.5 

758.0 

-1.58 

3 

5.24 

6.90 

1.50 

0.2 

0.9220 

22.1 

1.5 

758.0 

-1.66 

4 

8.00 

9.62 

1.60 

0.4 

1.2691 

22.1 

1.5 

758.0 

-1.62 

5 

11.02 

12.71 

1.70 

0.4 

0.9375 

22.1 

1.5 

758.0 

-1.69 

6 

12.81 

14.58 

1.80 

0.6 

1.2150 

22.1 

1.5 

757.9 

-1.77 

7 

14.51 

16.30 

1.90 

0.6 

1.0770 

22.1 

1.5 

757.9 

-1.79 

8 

17.58 

19.36 

2.00 

0.8 

1.1919 

22.1 

1.5 

757.9 

-1.78 

9 

19.58 

21.39 

2.10 

0.8 

1.0713 

22.1 

1.5 

757.9 

-1.81 

10 

22.72 

24.57 

2.20 

1.0 

1.1540 

22.0 

1.5 

757.9 

-1.85 

11 

25.78 

27.47 

2.30 

1.0 

1.0190 

22.0 

1.5 

758.0 

-1.69 

12 

6.11 

7.81 

63.00 

0.4 

1.0188 

22.1 

1.5 

758.0 

-1.70 

13 

11.65 

13.47 

67.70 

0.6 

0.8207 

22.1 

1.5 

758.0 

-1.82 

NO.  DATA  SETS  IS  13 

COEF’S  FOR  'DELTA  P=A+B*FLOW-»-C*  (FLOW*FLOW)  ' 

OR  'FLOW*(-B+SQR<Bt2-4*C* (A-DELTA  P)>)/(2*C)' 

WHERE  DELTA  P IS  IN  MILLIBARS,  FLOW  IS  IN  C.C./SEC 


ZERO  COEF  * A *-6.9785E-03 
FIRST  DEG  COEF  = B * 3.3637E+00 
SECOND  DEG  COEF=  C =-3.0688E-03 


DELP(DATA) 

FLOW (DATA) 

FLOW (CALC) 

DEVIATION 

% DVATION 

O.OOOOE'i'OO 

1.6666E-08 

2.0746E-03 

-2.0746E-03 

-1.2447E+07 

7.9506E+00 

2.3726E+00 

2.3708E+00 

1.8592E-03 

7.8360E-02 

1.1828E+01 

3.5548E+00 

3.5298E+00 

2.4930E-02 

7.0132E-01 

1.7168E+01 

5.1467E+00 

5.1298E-I-00 

1.6833E-02 

3.2708E-01 

2.3121H-^01 

6.9403E+00 

6.9194E+00 

2.0892E-02 

3.0102E-01 

2.6687E+01 

8.0111E+00 

7.9942E+00 

1.6919E-02 

2.1120E-01 

3.0020E+01 

9.0144E+00 

9.0005E+00 

1.3930E-02 

1.5453E-01 

3.5993E+01 

1.0819E+01 

1.0808E+01 

1.0133E-02 

9.3659E-02 

3.9919E+01 

1.2002E+01 

1.2001E+01 

1.7258E-03 

1.4378E-02 

4.6090E+01 

1.3874E+01 

1.3879E-I-01 

-5.6695E-03 

-4.0864E-02 

5.1899E+01 

1.5654E+01 

1.5654E+01 

-5.4929E-04 

-3.5089E-03 

1.3573E+01 

4.0246E+00 

4.0522E+00 

-2.7671E-02 

-6.8756E-01 

2.4496E+01 

7.2623E+00 

7.3335E+00 

-7.1266E-02 

-9.8132E-01 

SUM  OF  SQRS  OF  DEVIATION*  7.8122E-03 


DELP(CALC) 

FLOW (CALC) 

O.OOOOE+00 

2.0746E-03 

5.0000E-t-00 

1.4905E+00 

l.OOOOE-^01 

2.9830E+00 

1.5000E+01 

4.4796E+00 

2.0000E-I-01 

5.9804E+00 

2.5000E+01 

7.4853E+00 

3.0000E+01 

8.9944E-»'00 

3.5000E+01 

1.0507E+01 

4.0000E+01 

1.2025E+01 

4.5000E+01 

1.3547E+01 

5.0000E+01 

1.5073E+01 

Figure  2. 

Program  output  from  thermal  printer  for 

FULL  RANGE 
.4895E-02 
.2310E-02 
.9916E-01 
.0200E-01 
.5070E-01 
.0303E-01 
.6716E-01 
.2159E-01 
.0710E-02 
.8034E-02 
.5915E-03 
.3206E-01 
.5520E-01 


direction. 
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Figure  3.  Program  output  from  plotter. 


4.  USE  OF  PROGRAM  OUTPUT 

The  flow  meter  calibration  coefficients  calculated  by  this  prograun 
provide  a relationship  between  the  volume  flow  rate  at  the  midpoint  of 
the  meter  and  the  measured  pressure  drop  across  its  signal  ports.  The 
relationship  can  be  expressed  as  equation  (12a) : 


APg  = A + BQ^  + , 

by  which  it  is  assumed  that  the  gas  is  air  and  that  it  is  at  stamdard 
conditions  (20^C,  1 bar)  at  the  midpoint  of  the  meter.  If  steuidard 
conditions  do  not  exist  at  the  meter  midpoint,  then  APg  must  be 
calculated  by  use  of  equation  (11a)  before  cam  be  found  from  either 
equation  (12a)  or  from  the  program's  plot  output. 

If  the  flow  of  a gas  other  them  air  is  measured,  the  values  of  B amd 
C must  be  modified  to  reflect  the  resulting  chamges  in  y and  p , since 
from  equation  (12)  , 


B = by  and  C * cp 
s s 

Therefore,  eqiiation  (12a)  is  modified  to  become 

AP  = A + BiQ  + CiQ2  , (14) 

s m m 

where 

Bi  = B • — ^ = B * 

^s 

AIR 

y is  the  viscosity  of  the  gas  being  used  (as  given  in  teible  I) , 
s 

g 

C = C • (S.G.)  , 

S.G.  is  the  specific  gravity  of  the  gas  being  used  (as  given  in 
table  I) . 

Finally,  if  the  equivalent  volume  flow  at  conditions  other  than 
those  that  exist  at  the  midpoint  of  the  meter  (e.g.,  standard 
conditions)  is  desired,  then  equation  (13)  can  be  used. 


181.5  X 10~'7 
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APPENDIX  A. — PROGRAM  LISTING 


60  PRINT 
70  PRINT 


80  PRINT 


10  REM  *****CALFL0**2O  MAR  75**M.FUMKE********************* 

16  RF.M  *****IJSE  TO  CALIBRATE  VOL-O-FLO  METERS************** 

20  REM  *****E0  FOR  FIT  IS  LFAST  SOHARES  AND  SECOND  ORDER*** 

30  DIM  P(31 ).0(31 ).L(31 ),R(31 ),A(2,3),U(3).T(2).B(2.3) ,C(3).K(3) 
40  SELECT  PRINT  005(64) 

50  PRINT  "INSTRUCTIONS:" 

1.  MAX  NO.  DATA  SETS  ENTERED  IS  30." 

2.  DATA  SETS  SHOULD  BEGIN  WITH  'AB'  DIRECTION  DATA, 
TERMINATED  BY  'LEFT'  = 9.999,  FOLLOWED  BY  'BA' 
DATA,  TERMINATED  BY  'LEFT'  * 9.999." 

3.  DATA  CORRECTIONS  ARE  ALLOWED  AFTER  ALL  DATA  FOR 
CURRENT  DIRECTION  HAVE  BEEN  ENTERED." 

90  PRINT  " 4.  ON  PLOTTER,  SET  ZERO  REFERENCE  IN  LOWER  LEFT 
CORNER;  SET  SCALE  AT  X=Y=+8  INCHES." 

100  STOP  " KEY  CONTINUE,  EXECUTE" 

no  PRINT  HEX(03) 

120  INPUT  "DATE  (DAY  MONTH  YR) 

130  INPUT  "METER  S/N 
140  INPUT  "METER  RANGE(LPM) 

R*1 000/60 
150  D1S*"AB" 

160  SELECT  PRINT  015(80) 

170  PRINT  :PRINT 
180  PRINT  "DATE  (DAY  MONTH  YR) 

190  PRINT  "METER  S/N  ";  S 

200  PRINT  "METER  RANGE 
ITHEN  2A0 
lOTHEN  250 

lOOTHEN  260:G0SUB  300:G0T0  310 
240  ROSUB  270: GOTO  310 
250  GOSUB  280:G0T0  310 
260  GOSUB  290: GOTO  310 
270  PRINTUSING  1 840, R; : RETURN 
280  PRINTUSING  1860,R; :RETURN 
290  PRINTUSING  1880,R; :RETURN 
300  PRINTUSING  1 900, R;: RETURN 
310  PRINT  " CCPS  (";R$;: PRINT  " LPM)" 

320  REM  *****NORMALIZE  AND  PRINT  DATA*** 


D$ 

S 

R$: CONVERT  R$TO  R:R* 


ns 


210  IF  R 
220  IF  R 
230  IF  R 


330  FOR  K=1T0  3:K(K)*0:NEXT  K 

340  L-0:  M»2:  PRINT  :PRINT  :PRINT  :PRINT  "FLOW  DIRECTION  IS 
D1$:PRINT 

350  PRINTUSING  360, "NO" , "LEFT" , "RGHT" , "P2" ,"VOL" , "TIME" , "T2" 
,"PO","BARP","DIFF" 

360  X ##  §m  mp  ##  ###  ####  MPP  ##  #### 
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370  L(l),R(l),B,Tl,Pn=0:V=1E-n:I,T=l:Bl=760:I=l:N3=n 
3R0  D=L(I)-R(I):PRINTIISING  400,1, L(I) ,R(I),B,V,T,T1 ,P0,B1 ,D 
300  SELECT  PRINT  005:PRINT  HEX(03) :PRINT  "DATA  POINT" ; I+l : SELECT 
PRINT  015(80) 

400  % PH  -PP.PP  -»pj  -P.mp  -pp.p  -#.#  -###.# 

410  P{I)=(L(I)+R(I))*2.488V?.5A 

420  0( I )=1 000*V/T/60* (1.871 7*P0+B1 )/ ( .7521 8*P ( I )/2+7 . 521 8*B+B1 ) 
430  P(I)=P(I)/(1+9/5*.00146*(T1-20)+.0013/760*(752.18-.75218*P(I 
)/2-7.5218*B-B1)): 1=1+1: IF  N2=1THEN  550:N9=I 
440  INPUT  "LEFT  MANOMETER  READING(CM.W.C.)=",L(I) 

450  IF  L(I)=9.999THEN  550 

460  INPUT  "RIGHT  MANOMETER  READING(CM.W.C.)=",R(I) 

470  INPUT  "BACK  PRESSURE (KPA)=",B 

480  INPUT  "VOLUME(LITERS)=",V 

490  INPUT  "TIME(MIN)=",T 

500  INPUT  "TEMPERATURE (DEG  C)=",T1 

510  INPUT  "VOLUME  PRESSURE(IN.W.C.)=",PO 

520  INPUT  "BAROMETRIC  PRESSURE(MM.HG. )=",R1 -.GOTO  380 

530  REM  *****DATA  CORRECTIONS*** 

540  N8=0 

550  N2=0: INPUT  "DO  YOU  WANT  TO  MAKE  DATA  C0RRECTI0NS(1=YES)", 

N2 

560  IF  N2=1THEN  570; IF  N3=1THEN  880:GOT0  1120 

570  IF  N8=1THEN  580:PRINT  :PRINT  "POINTS  CORRECTED  INCLUDE" 

580  N8=1: INPUT  "WHICH  DATA  POINT  NUMBER", I: GOTO  440 
590  REM  *****PRINT  COEFF'S*** 

600  SELECT  PRINT  015(80) 

610  PRINT  rPRINT 
620  Z=0:F0R  K=1T0  3 

630  IF  K(K)riTHEN  640;IF  K(K)H-1THEN  640:IF  K(K)r]INT(K(K))THEN 
640: GOTO  650 
640  K(K)=0 

650  IF  K(K)=OTHEN  670; IF  Z=1THEN  660:PRINT  "POINTS  DELETED  INCLU 
DE":Z=1 

660  PRINT  K(K) 

670  NEXT  K 

680  PRINT  "NO.  DATA  SETS  IS  ";  I-1-Z1 

690  PRINT  :PRINT  " COEF'S  FOR  'DELTA  P=A+B*FLOW+C*(FLOW*FLOW) 
:PRINT  " OR  'FL0W=(-B+S0RrB!2-4*C*(A-DELTA  P)])/(2* 

O'" 

700  PRINT  :PRINT  " WHERE  DELTA  P IS  IN  MILLIBARS,";: 

PRINT  " FLOW  IS  IN  C.C./SEC" 

710  PRINT 

720  PRINT  "ZERO  COEF  = A =";:PRINTUSING  730,  C(1):PRINT 
"FIRST  DEG  COEF  = B =";;PRINTUSING  730,  C(2):  PRINT 
"SECOND  DEG  COEF=  C =";;  PRINTUSING  730,  C(3) 
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730  «-#.####!!!! 

740  REM  *****PRINT  DATA  ANALYSIS*** 

750  PRINT  tPRINTUSINR  760,"DELP(0ATA)"."FL0W(DATA)",''FL0W(CALC)" 
."DEVIATION","?  OVATION","?  FULL  RANGE" 

760  ? mnmm  mmm  ###### 

770  S1,Z=0:F0R  N=1T0  I-1:F0R  K=1T0  3: IF  K(K)riNTHEN  780:2=1 
780  NEXT  K:IF  Z'lTHEN  830 

790  X1=(-C(2)+S0R(C(2)*C(2)-4*C(3)*(C(1)-P(N))))/(2*C(3)) 

800  W=0(N)-X1:  S1=S1+H*W:  W0=100*W/Q(N):  W5=100*W/R 
810  PRINTUSING  820,  P(N),0(N) ,X1 ,W,W0,W5 
820  ?-#.^###!!!!  -#.^##«!!!!  -#.#### 
!!!!  -#.#W!!!! 

830  Z=0:NEXT  N:  PRINT 

840  PRINT  " SUM  OF  SQRS  OF  DEVIATION="; : PRINTUSING  850,  SI 

850  ?-#.##4#mm 

860  N3=1:G0T0  530 

870  REM  *****DATA  DELETIONS*** 

890  N4=0: INPUT  "DO  YOU  WANT  TO  DELETE  A DATA  P0INT(1=YES)".N4 
900  IF  N4=1THEN  910:IF  N3=2THEN  1120:IF  N8=1THEN  1120:G0T0  930 
910  K=K+1:N3=2 

920  INPUT  "WHICH  DATA  POINT  NUMBER", K(K) : IF  Kr3THEN  890:G0T0  112 
0 

930  IF  D1S="BA"THEN  970 

940  REM  *****SET  PRESSURE  AND  FLOW  SCALE*** 

950  04=2*R 

960  FI =625/50: F2=625/04 

970  REM  *****GENERATE  PRESS-FLOW  TABLE*** 

980  PRINT  : PRINTUSING  990,  "DELP(CALC)","FLOW(CALC)" 

990  % 

1000  IF  D1$="BA"THEN  1010:G0SUB  '18:G0T0  1020 

1010  GOSUB  '19 

1020  FOR  N=0  TO  10:  Y=5*N 

1030  X1=(-C(2)+SQR(C(2)*C(2)-4*C(3)*(C(1)-Y)))/(2*C(3)) 

1040  PRINTUSING  1050.  Y,X1 
1050  ?-#.#4##!!!!  -#.44##!!!! 

1060  REM  *****PLOT  PRESS-FLOW  CURVE*** 

1070  IF  NiOTHEN  1080:G0SUB  '1 1 (Y,X1 ):GOTO  1090 
1080  GOSUB  '12(Y,X1) 

1090  NEXT  N 

1100  SELECT  PRINT  005 (64): GOTO  1280 
1110  REM  *****COEFF  CALCULATION*** 

1120  10,11 ,I2,K0,K1 ,K2,K3,K4,K5,K6,K7=0:W=1 : I*N9:Z1 ,Z=0 
1130  FOR  N=1T0  I-1:F0R  K=1T0  3:IF  K(K)r]NTHEN  1140:Z=1 :Z1=Z1+1 
1140  NEXT  K:IF  Z=1THEN  1160 
1150  IO=IO+W*0(N):KO=KO+P(N)*W:K4=K4+W 
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1160  Z=0:NEXT  N:K0=K0/K4:I0=I0/K4 

1170  FOR  N=1T0  I-1:F0R  K=1T0  3:IF  K(K)r]NTHEN  11R0:Z=1 

1180  NEXT  K:IF  Z=1THEM  1200 

1190  H1=0(N)-IO:I1=I1+0(N)*H1*H1*W:K1=K1+P{N)*H1*W:K5=K5+H1*H1*W 

1200  Z=0:NEXT  N:K1=K1/K5:J1=K5/K4:I1=I1/K5 

1210  FOR  N=1T0  I-1;F0R  K=1T0  3:IF  K(K)nNTHEN  1220:Z*1 

1220  NEXT  K:IF  Z=1THEN  1240 

1230  H2=(Q(N)-I0)*(Q(N)-lf )-Jl : I2=I2+Q(N)*H2*H2*W:K2=K2+P{N)*H2 
*W: K6=K6+H2*H2*W 

1240  Z=0:NEXT  N:K2*K2/K6:J2=K6/K5:I2=I2/K6 
1250  C(1)=K0-K1*I0+K2*{I0*I1-J1) 

1260  C(2)=Kl-K2*(I0+Il):C(3)=K2:fi0T0  590 
1270  REM  *****SUBROUTINE  FOR  PLOTTING*** 

1280  V3=5:  V4=04/5 
1290  GOSUB  '20(V3,V4) 

1300  V4=50:  V5=10:  V7=04:  V8=04/5 

1310  GOSUB  '21(V4,V5,V7,V8):PL0T  [,,R]:GOTO  2000 

1320  REM  *****ROUNDOFF  COMPENSATION*** 

1 330  DEFFN ' 1 0 ( X9 . Y9 ) : X8= INT( FI *X9-X0+ . 5 ) : Y8= I NT { F2*Y9-Y0+ . 5 ) : 
X0=X0+X8 : Y0=  Y0+Y8 : RETURN 
1340  REM  *****PLOT  POINT  UP*** 

1350  DEFFN' 11(X9,Y9): GOSUB  '10(X9,Y9):PL0T  rX8,Y8,U]:RETURN 
1360  REM  *****PLOT  POINT  DOWN*** 

1370  DEFFN'12(X9,Y9):G0SUB  '10(X9,Y9):PLOT  rX8,Y8,  Pi -.RETURN 
1380  REM  *****PLOT  SYMBOL*** 

1390  DEFFN '17{X9,Y9):G0SUB  '10(X9,Y9):PL0T  rX8,Y8,F$]:RETURN 
1400  REM  *****GOTO  ORIGIN  OF  'AB'  PLOT*** 

1410  DEFFN '1 8: X0,Y0=0: PLOT  r.,R]. [125,125,0]:  RETURN 
1420  REM  *****GOTO  ORIGIN  OF  'BA'  PLOT*** 

1430  DEFFN '1 9: X0,Y0=0: PLOT  [,,Rl, [500,62, U],[500,63,U]:  RETURN 
1440  REM  *****DRAW  AXIS*** 

1450  DEFFN'20(V3,V4) 

1460  IF  D1S*"AB"THEN  1480 
1470  GOSUB  '19-.G0T0  1490 
1480  GOSUB  '18 
1490  V5=0 

1500  GOSUB  '12(V5,0) 

1510  PLOT  [,10,D],[,-20,D],[,10,D] 

1520  V5=V5+V3 

1530  IF  V5[=50  THEN  1500 

1540  V5=0 

1550  GOSUB  ’12(0,V5) 

1560  PLOT  [10,,D],[-20,,D],[10,,D] 

1570  V5=V5+V4 

1580  IF  V5[»04  THEN  1550 

1590  RETURN 

1600  REM  *****SCALE  AXIS*** 
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1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

0. 


DEFFN'21(V4,V5,V7,V8) 

IF  D1$«"AB''THEN  1640 
GOSIJB  '19:G0T0  1650 
GOSUB  *18 
SELECT  PRINT  413 
GOSUB  1940 
V9-0 

GOSUB  '10(V9,0) 

X8-X8-39:  X0«X(H28;  Y8*Y8-20:  Y0=Y0-20 
PLOT  [X8,Y8,U] 

0*V9:  GOSUB  1970 
V9»V9+V5 

IF  V9[=V4  THEN  1680:  IF  D1$-"AB"THEN  1740:  GOSUB 
1750 

GOSUB  *18 
V9=0 

PLOT  [-120,0.U1 

0»V9:IF  0«0THEN  1810:IF  0[1THEN  1800 
IF  OriOTHEN  1810 

IF  OClOOTHEN  1820:G0SUB  1890:G0T0  1910 

GOSUB  1830:G0T0  1910 

GOSUB  1850: GOTO  1910 

GOSUB  1870: GOTO  1910 

PRINTUSING  1840,0;: RETURN 

%-,m 

PRINTUSING  1860,0; :RETURN 

PRINTUSING  1880,0;  :RETIIRN 

PRINTUSING  1900,0;: RETURN 

PLOT  r-110,125,U]:V9»V9+V8 
IF  V9[-V7  THEN  1770 
SELECT  PRINT  005:RETURN 
PLOT  [2,,C],[22,,S]:  RETURN 
PRINTUSING  1960,0;:  RETURN 

%-m 


'19:G0T0 


PRINTUSING  1980,0; :RETURN 
%-## 

REM  *****LABEL  AXIS*** 

PLOT  [2,,C],[25,,S],[,,R] 

SELECT  PRINT  413 

PLOT  [225,40,U]:IF  D1$*"AB''THEN  2030: PLOT  2[438,0,U] 

PRINT  "DELTA  P (MILLIBARS)" 

PLOT  [,,Rl,ri90,775,U]:IF  D1$*"AB"THEN  2050:PL0T  2[438,0,U] 
PRINT  "FLOW  (CC/SEC)":IF  D1$="AB"THEN  206O:PL0T  [390,361 
,U]:PRINT  ";";I-1-Z1;:PL0T  [-20,0,U]:PRINT  "'BA'":PLOT  [-85 
-185,II]:PRINT  " 'BA‘  DIRECTION": GOTO  2260 
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2060  R6M  *****LABEL  PLOT*** 

2070  PLOT  [-500,0,Ul,r60,460,Ul 
2080  PRINT  "METER  S/N" 

2090  PLOT  r-225,-45,Ul 
2100  PRINT  "METER  RANGE" 

2110  PLOT  [-275,-45,01 

2120  PRINT  "NO.  OF  DATA  POINTS" 

2130  PLOT  [-450,-45,0] 

2140  PRINT  "DATE" 

2150  PLOT  [450,-48,0] 

2160  PRINT  "COEFF'S  FOR  P=A+BQ+CO"; :PL0T  [-35,10,0],[1 ,,C]:PRINT 
2:PL0T  [-650,-50,0], [2, ,C]:PRINT  "OR  Q=(-B+SQR(B"; :PLOT 
[-30,10,0],[l,,Cl:PRiNT  2:PL0T  [-30,-10,0],[2,,C]:PRINT 
"-4C(A-P)))/2C" 

2170  PLOT  [275,225,0]:PRINT  S:PL0T  [-150,-45,0] 

2180  IF  R[1THEN  2210 
2190  IF  R[10THEN  2220 

2200  IF  R[100THEN  2230:G0S0B  300:G0T0  2240 
2210  GOSOB  27O:G0T0  2240 
2220  GOSOB  280: ROTO  2240 
2230  GOSOB  290 

2240  PRINT  " CCPS":PLOT  [-250,-45,0]:PRINT  I-1-Z1;:PL0T  [-20,0,0 
]:  PRINT  "'AB"':PLOT  [-160,-45,0]:PRINT  D$:PLOT  3[-505,-46,U] 

2250  PRINT  " 'AB'  DIRECTION" 

2260  PLOT  [-275,-50,01 

2270  PRINT  "A  =";:PRlifrOSING  730,C(1) 

2280  PLOT  2[-175,-20,0] 

2290  PLOT  r2,,Cl,[25,,S] 

2300  PRINT  "B  tPRINTOSINR  730, C(2) 

2310  PLOT  2[-175,-20,0] 

2320  PLOT  [2,,C],[25,,S] 

2330  PRINT  "C  ="; rPRIfTTOSING  730,C(3):PL0T  [,,R]:SELECT  PRINT 
005 

2340  F$="+":PL0T  [1 ,,C],[0,,S]:IF  D1$="BA"THEN  2350:G0S0B  '18: 
GOTO  2360 
2350  GOSOB  '19 

2360  FOR  N=1T0  I-1:F0R  K=1T0  3:IF  K(K)[]NTHEN  2370:Z=1 
2370  NEXT  K:IF  Z=1THEN  2390 
2380  GOSOB  '17(P(N),Q(N)) 

2390  Z=0:NEXT  N:IF  DlS="BA"THEN  2400 :D1$=" BA": SELECT  PRINT  015( 
80): GOTO  320 

2400  SELECT  PRINT  005: PLOT  [,,R]:EN0 
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DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION,  BUILDING  5 
ALEXANDRIA,  VA  22314 
ATTN  DDC-TCA  (12  COPIES) 

OFC,  CHIEF  OF  RESEARCH  & DEVELOPMENT 
USA  RSCH  & DEV  GROUP  (EUROPE) 

BOX  15 

FPR  NEW  YORK  09510 
ATTN  LTC  EDWARD  E.  CHICK 

CHIEF,  MATERIALS  BRANCH 

COMMANDER 

US  ARMY  MATERIEL  DEVELOPMENT 
& READINESS  COMMAND 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
ATTN  DRCRD,  RES,  DEV, 

& ENGR  DIRECTORATE 

COMMANDER 

USA  ARMAMENT  COMMAND 
ROCK  ISLAND,  IL  61201 
ATTN  DRSAR-ASF,  FUZE  DIV 
ATTN  DRSAR-RDF,  SYS  DEV  DIV  - FUZES 

COMMANDER 

USA  MISSILE  & MUNITIONS 
CENTER  & SCHOOL 
REDSTONE  ARSENAL,  AL  35809 
ATTN  ATSK-CTD-F 

COMMANDER 
PICATINNY  ARSENAL 
NUCLEAR  DEVELOPMENT  & ENGR.  DIR 
DOVER,  NJ  07801 
ATTN  SARPA,  MR.  A.  FRANZ 
ATTN  SARPA-ND-C-C , MR.  DAVID  SAMPAR 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
WASHINGTON,  DC  07801 
ATTN  AIR  52022A1, 

MR.  WILLIAM  LAUGHREY 

NASA,  LANGLEY  STATION 
HAMPTON,  VA  23365 
ATTN  MS  494,  MR.  H.  DOUGLAS  GARNER 

COMMANDER 

AIR  FORCE  FLIGHT  DYNAMICS  LAB 
FLIGHT  CONTROL  DIVISION  (FDCL) 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433 
ATTN  MR.  HARRY  M.  SNOWBALL 

COMMANDER 

US  NAVAL  AIR  DEVELOPMENT  CENTER 
AEROMECHANICS  DEPARTMENT 
WARMINSTER,  PA  18974 
ATTN  MR.  HORACE  WELK 


COMMANDER 

US  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 
ATTN  DRSMI-RGC,  MR.  JAMES  WILLIAMS 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  DC  20360 
ATTN  CODE  AIR  52022A,  MR.  JOHN  BURNS 

C0MM7VNDER 

NAV7VL  SHIP  SYSTEMS  COMMAND 
WASHINGTON,  DC  20360 
ATTN  SHIPS  03414,  MR.  ARTHUR  CHAIKIN 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 
WRIGHT-PATTERSON  AIR  FORCE  BhSE,  OH  45433 
ATTN  AFIT/ENM,  MR.  M.  E.  FRANKE 

COMMANDER 

US  ARMY  WEAPONS  COMMAND 
ARMORED  WEAPONS  SYSTEMS  DIR 
ROCK  ISLAND,  IL  61201 
ATTN  SWERR-T-AC,  LORRAINE  D.  WRIGHT 

COMMANDER 

USAMERDC 

FT.  BELVOIR,  VA  22060 
ATTN  SMEFB-EP,  MR.  ROBERT  WARE 

HARRY  DIAMOND  LABORATORIES 
ATTN  MCGREGOR,  THOMAS,  COL,  COMMANDING 
OFFICER/FLYER,  I. N. /LANDIS,  P.E./ 
SOMMER,  H. /CONRAD,  E.E. 

ATTN  CARTER,  W.W.,  DR.,  ACTING  TECHNICAL 
DIRECTOR/MARCUS,  S.M. 

ATTN  KIMMEL,  S.,  10 

ATTN  CHIEF,  0021 

ATTN  CHIEF,  0022 

ATTN  CHIEF,  LAB  100 

ATTN  CHIEF,  LAB  200 

ATTN  CHIEF,  LAB  300 

ATTN  CHIEF,  LAB  400 

ATTN  CHIEF,  LAB  500 

ATTN  CHIEF,  LAB  600 

ATTN  CHIEF,  DIV  700 

ATTN  CHIEF,  DIV  800 

ATTN  CHIEF,  LAB  900 

ATTN  CHIEF,  LAB  1000 

ATTN  RECORD  COPY,  BR  041 

ATTN  HDL  LIBRARY  (3  COPIES) 

ATTN  CHAIRMAN,  EDITORIAL  COMMITTEE 

ATTN  CHIEF,  047 

ATTN  TECH  REPORTS,  013 

ATTN  PATENT  LAW  BRANCH,  071 

ATTN  MCLAUGHLIN,  P.W. , 741 

ATTN  CHIEF,  310  (15  COPIES) 

ATTN  CHIEF,  340  (10  COPIES) 

ATTN  CAMPAGNUOLO,  C.,  940 
ATTN  FUNKE,  M.  F.  (5  COPIES) 
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